Nuclear prolate-shape dominance with the Woods-Saxon potential 
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We study the prolate-shape predominance of the nuclear ground-state deformation by calculating 
the masses of more than two thousand even-even nuclei using the Strutinsky method, modified 
by Kruppa. and improved by us. The influences of the surface thickness of the single-particle 
potentials, the strength of the spin-orbit potential, and the pairing correlations are investigated 
by varying the parameters of the Woods-Saxon potential and the pairing interaction. The strong 
interference between the effects of the surface thickness and the spin-orbit potential is confirmed to 
persist for six sets of the Woods-Saxon potential parameters. The observed behavior of the ratios 
of prolate, oblate, and spherical nuclei versus potential parameters are rather different in different 
mass regions. It is also found that the ratio of spherical nuclei increases for weakly bound unstable 
nuclei. Differences of the results from the calculations with the Nilsson potential are described in 
detail. 

PACS numbers: 21.10.Ft, 21.10.Gv, 21.10.Ky, 21.10.Pc 



I. INTRODUCTION 

It has long been known that most of the deformed 
nuclei have prolate rather than oblate shapes [l|. This 
prolate-shape dominance is the subject of this paper, in 
which we examine the effects of several important factors 
to increase or decrease such preference for prolate shapes. 

The properties of a nucleus are completely determined 
by the Hamiltonian for quantum mechanical many-body 
systems consisting of nucleons. From this point of view, a 
possible direction of research for the origin of the prolate- 
shape dominance is to change artificially each term of 
the Hamiltonian and examine its effect on the prolate- 
shape dominance. The most successful outcome of such 
a study would be a discovery of a direct correspondence 
between the prolate-shape dominance and a specific term 
of the elementary nucleon-nucleon interaction, in an anal- 
ogous fashion as the tensor force causes a mixture of the 
d wave in the wavefunction of a deuteron. However, no- 
body knows whether there exists such a simple relation. 

In mean-field approximations all the influences of the 
many-body Hamiltonian on the nuclear shape are con- 
veyed through the mean-field (one-body) Hamiltonian. 
Hence, one can seek the origin of the prolate-shape dom- 
inance in the properties of the mean field, instead. If 
one also adopts the macroscopic-microscopic method Q , 
which uses only the single-particle level density, one can 
determine the nuclear shape starting from the mean-field 
Hamiltonian without referring to the underlying many- 
body Hamiltonian. This is the choice of this paper. One 
can reasonably expect that the linkage from the mean- 
field Hamiltonian to the nuclear shape may be simpler 
than those from the many-body Hamiltonian. Moreover, 



the former linkage has an established route to understand 
in terms of classical periodic orbits 0, 3 , while the latter 
has none presently. 

We have been seeking the origin of the prolate-shape 
dominance in the one-body potentials in the macroscopic- 
microscopic framework. In Ref. [f| , we employed the Nils- 
son potential @ and calculated the ratio of prolate-shape 
nuclei among ~1800 even-even nuclei as a function of 
the strengths of the I 2 term and the spin-orbit term. We 
found that the ratio of prolate nuclei oscillates versus the 
strength of the spin-orbit term when the I 2 term has the 
standard strength. This oscillation disappears when the 
I 2 potential is weaker. The amplitude of the oscillation 
is quite large. The ratio of prolate-shape nuclei oscillates 
between ^40% (at ±40% strengths of the standard) and 
>80% (at ±100% and 0% strengths). This result indi- 
cates a strong interference between the effects of the two 
terms of the potential. 

However, this interference might be exaggerated be- 
cause of the affinity of the two terms in the sense that 
both of them include the orbital angular momentum op- 
erator I. The I 2 term is an approximation to the radial 
profile of the central potential like a square well. In a let- 
ter , we replaced the Nilsson potential with the Woods- 
Saxon (WS) potential to eliminate the I 2 potential. As a 
result, we still find the oscillation, although its amplitude 
was decreased by a factor ~ i. 

Employment of the WS potential required various 
modifications of the method Q. First, the WS poten- 
tial has a continuum spectrum, which cannot be treated 
by the standard Strutinsky method. This problem was 
solved by the Kruppa's prescription Q to subtract the 
level density of the free-particle spectrum. Concerning 
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the remaining problem of the plateau condition, we in- 
vented the reference density method [§[ to make the shell 
energy insensitive to the smearing width. Second, the 
BCS equation with a constant pair-scattering matrix el- 
ement G also suffers from the continuum spectrum. Wc 
extended the Kruppa's prescription to the paring to sub- 
tract the contributions from free-particle spectrum in the 
BCS equation. This treatment also improves the calcu- 
lation of the nuclear radius fioj . Third, available param- 
eter sets of the WS potentials are not consistent with 
the macroscopic part especially for nuclei far from the 
/3-stability line. We presented a method based on the 
Thomas-Fcrmi approximation to adjust the depth of the 
central potential to reproduce the drip lines predicted by 
the macroscopic energy. 

In this paper, we show and discuss the results of the 
calculations with the WS potential in much more detail 
than in the letter . We try not a single but six parame- 
ter sets of the WS potential. We confirm the convergence 
of the results versus the size of the truncated oscillator 
basis for the diagonalization of the single-particle Hamil- 
tonian. We show that different definitions of the ratio of 
prolate-shape nuclei do not change the conclusions. We 
also consider the effects of changing the strength of the 
pairing interaction, which we studied so far only for the 
Nilsson potential in Ref. We show the distribution 
of prolate and oblate nuclei in the nuclear chart for vari- 
ous modified potentials and confirm the independence of 
the main conclusions from the choice of potentials. 

In the rest of the introduction, we survey the earlier 
works. 

For the anisotropic harmonic oscillator potential, 
which is the most basic approximation to the one-nucleon 
potential, prolate (oblate) shapes tend to have lower en- 
ergies than oblate (prolate) ones when the major shell is 
less (more) than half filled because the particles (holes) 
in the low Vt orbitals drives strongly the nuclear shape 
toward prolate (oblate) direction. Here, fi is the projec- 
tion of the single-particle angular momentum onto the 
symmetry axis of the mean-field potential. This suggests 
an approximately equal number of prolate and oblate nu- 
clei Indeed, our quantitative estimation Q showed 
that the ratio of prolate nuclei among deformed nuclei is 
55% for a harmonic oscillator potential, while the exper- 
imental ratio of prolate nuclei seems much larger. 

The earliest attempt to explain the prolate-shape pref- 
erence is a work by Lemmer and Weisskopf [l2j|. A term 
oc r 4 was added to the harmonic oscillator potential to 
steepen the wall. Then it overrode the normal tendency 
of the oscillator potential to deform into an oblate shape 
for a more-than-half-filled major shell. It suggests that 
the origin of the prolate-shape dominance is the radial 
profile of the potential. 

However, there was also an opinion that the origin is 
the spin-orbit potential. For example, Ref. [l[ states that 
it is the spin-orbit potential which breaks the even sit- 
uation by weakening the oblate-shape shell effect in sd 
shell nuclei. 



The truth has turned out [f| that both of the two 
salient features of the nuclear single-particle potential, 
the square-well like radial profile and the spin-orbit po- 
tential, play essential roles to give rise to the prolate- 
shape dominance. The dominance can be equally repro- 
duced with and without the spin-orbit potential but it 
does not deny the participation of the spin-orbit poten- 
tial, because half-strength spin-orbit potential destroys 
the dominance. 

The relation between potential and spectrum can be 
understood in terms of classical periodic orbits. Putting 
aside the spin-orbit potential, Frisk analyzed classical pe- 
riodic orbits in an ellipsoidal cavity and found that the 
radial dependence is an origin of the prolate-shape dom- 
inance [3| ■ The strength of the shell effect at the Fermi 
surface changes strongly in the prolate side while it stays 
almost constant in the oblate side as a function of the 
magnitude of deformation. As the classical periodic or- 
bit responsible for this asymmetry between prolate and 
oblate shapes, Frisk payed attention to a triangular orbit 
in the meridian plane whose period changes little in the 
oblate side owing to the volume conservation condition. 

Arita used a power-law potential oc r a Q , which can 
interpolate continuously between the square-well poten- 
tial and the harmonic oscillator potentials, for the same 
kind of analysis and obtained essentially the same con- 
clusion as Frisk [|| . Concerning the classical periodic or- 
bits responsible for the prolate-shape dominance, Arita 
insisted the importance of bridge-orbit bifurcations. He 
is also planning to treat the spin-orbit potential in the 
periodic orbit theory fl3j . 

Hamamoto and Mottelson presented a plain explana- 
tion to the fact that the level density in the Fermi level 
changes little in the oblate side flij . In the Nilsson di- 
agram of single-particle spectrum, the degenerated lev- 
els belonging to a subshcll at spherical shape fans out 
(spreads) versus increasing magnitude of deformation in 
both prolate and oblate sides. However, this fanning out 
is small for levels with low A (the projection of orbital an- 
gular momentum onto the symmetry axis) in the oblate 
side because they are pushed down by many orbitals with 
the same quantum number belonging to upper major 
shells. Since the same spectrum leads to the same value 
of the shell effect, no energy can be gained by deforming 
into oblate shapes. 

However, Arita Q found a counterexample to this ex- 
planation. For a low-power potential oc r , the fanning 
is suppressed for low- A orbitals in prolate shapes, while he 
found no oblate-shape dominance. Therefore, the expla- 
nation in terms of fanning applies indeed to nuclcar-likc 
potentials but probably not to different types of poten- 
tials. We think that the reason for this limitation is that 
the explanation in terms of the fanning mechanism does 
not explicitly consider the volume conservation condition 
which played an essential role in Frisk's argument @ . 

Finally, we mention two other factors which have minor 
effects to prefer prolate shapes. 

One is the Coulomb interaction, which tends to clon- 
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gate the nuclear shape rather than flatten it in order to 
diminish positive electrostatic potential energy between 
protons. However, this effect does not seem to play a 
decisive role to give rise to the prolate-shape dominance 
because it increases as Z 2 and is strongest in heavy nu- 
clei while the actual dominance is already very clear in 
middle-weight nuclei. This effect will be revisited in our 
future paper for more detailed studies. 

The other is the angular-momentum projection of 
mean-field solutions into zero or low-spin spin states, 
which has a tendency to lower the energy of prolate 
states relative to that of oblate ones, because the for- 
mer have generally larger moments of inertia than the 
latter. This effect seems to have already been argued 
by Zickendraht [f5| as the difference of the volume ele- 
ment for the collective coordinates between prolate and 
oblate shapes, which is originated in the difference of the 
available configuration space in the spherical shell model. 
However, mean-field models and their phenomenologi- 
cal interactions are determined to reproduce the exper- 
iments without angular momentum projection and thus 
it is widely believed that they includes the effects of the 
projection. A quantitative examination of this belief will 
be a subject of our future study. 

In Sec. HU we explain about the potentials, our theo- 
retical framework, some practical information of our cal- 
culations, and definitions of the ratio of prolate nuclei. 
In Sec. IIII1 we give the results and discussions on them. 
In Sec. IIV1 we summarize the conclusions of this paper. 



II. METHODS 

We use the Woods-Saxon- Strutinsky [lj] and the 
Nilsson- Strutinsky methods @ to calculate the nuclear 
energy surface in the plane spanned by the axially sym- 
metric quadrupolc and hexadecapole deformation param- 
eters and to determine the shape of the ground state of 
even-even nuclei from Oxygen to superheavy elements 
between drip lines. For both methods, the single-particle 
Hamiltonian is a sum of a free kinetic energy and a po- 
tential energy, 

where m is the bare nucleon mass (the average of a pro- 
ton and a neutron masses). The Nilsson and the WS 
potentials used for V are described in Sees. Ill Al and III Bl 
respectively. The latter potential has a continuum part 
in the spectrum, which causes difficulties in the standard 
version of the Strutinsky method. Remedies we have cho- 
sen Q are summarized shortly in Sec. Ill Cl Precision and 
other practical aspects of numerical computations are de- 
scribed in Sees. Ill Dl and III El The quantity which we 
study intensively in this paper is introduced in Sec. Ill Fl 



A. The Nilsson potential 

The potential of the Nilsson model [3(| is expressed as 

V(r)=j {ulx 2 +u\y 2 +u 2 z 2 ) + h w° r 2 y^^o (*t) 
—IflsK^h Ldo h ■ S — fuKN/J-Nfi loq (If— (It 

(2) 

where uj± and lo z are related to a quadrupole deforma- 
tion parameter £2 through ui±=u>o (1 + ^2) and uj z ~uj^ 
(l — ^£2) with o?o determined by the condition of a vol- 
ume conservation uP'yUiz = u>q. Denoted by e 4 is a hexade- 
capole deformation parameter. Orbital and spin angular 
momenta are denoted as I and s, respectively and the 
subscript t means the usage of the stretched coordinates. 
The third term is a spin-orbit potential, while the fourth 
term is called the I 2 potential. The latter term propor- 
tional to the square of the orbital angular momentum is 
used to simulate roughly the change of the radial profile 
of the central potential from oc r 2 of the harmonic oscil- 
lator. With a standard negative strength of the fourth 
term (—fuKN/J-N) and without the spin-orbit potential, 
the spectrum is modified from that of the harmonic os- 
cillator toward that of the square- well potential. 

We adopt the standard values given in Table I of 
Ref. fl7j for the parameters kjv and fiN which are de- 
pendent on the harmonic oscillator shell N osc . By chang- 
ing the set of values of the multipliers (fi Sl fa) from the 
standard values (I, I), one can examine how the combi- 
nation of these two potentials affects the prolate-shape 
dominance. For example, // s =0 (1) means no (stan- 
dard) spin-orbit potential, while fu=Q (1) corresponds to 
harmonic-oscillator (standard) radial profile of the cen- 
tral potential. 

In this paper, we calculate the ratio of prolate nuclei 
for 31 x 16 = 496 kinds of artificially modified Nilsson 
potentials defined by the combination of values of the 
multipliers (fi s ,fu) in ranges —1.5 < fi s < 1.5 and < 
/„ < 1.5 with Af ls = Af u = 0.1. 

As a footnote, the eA-parametrization of nuclear de- 
formation in the Nilsson potential is different from the 
conventional /3\-paramctrization based on the expan- 
sion of nuclear radius in terms of the spherical harmon- 
ics [l8|. Including only the A = 2 deformation and re- 
quiring the ratio of axes is the same between the two 
kinds of parametrizations, as an example, one obtains 
€2 = 0.95 02 for small deformations while £2 — 1.14 /?2 for 
the ratio 2:1 of the ideal supcrdeformation. 

In this paper, we describe quantities related to shapes 
in term of /?2- For example, we give the thresholds to 
distinguish prolate or oblate shapes from the spherical 
shape in terms of the values of $2- For the solutions of 
the Nilsson-Strutinsky method, we apply these thresh- 
olds directory to £2, not calculating the value of /?2 cor- 
responding to the value of £2 , for the sake of simplicity. 
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B. The Woods-Saxon potential 

The Woods-Saxon (WS) potential is a finite-depth po- 
tential having a flat central and a steep surface parts. It 
is no doubt a better approximation to the central part of 
the nuclear potential V CE than the Nilsson potential. The 
spin-orbit potential V so is also reasonably expressed us- 
ing the gradient of this form of a potential. The Coulomb 
potential V co should be added for protons. Then, the 
sum of these three potentials, which we also call the WS 
potential, is expressed as 

V = V ce + V so + ^(1-t 3 )V co , (3) 

where T3 is the third component of the nucleon's isospin 
multiplied by 2 (1 for neutrons and —1 for protons), 



Vce — Vws( r , Voce, k cej Roce, /o"cei /3), 



fls^s 



2ra rc dC 



(4) 



(5) 



[VV' ws (r; V 0CE , k SO i Roso, /a 8 soi 

erx -\7 

i 



In Eq. ([5]), m rc d = ^-j^-m with m being the bare nucleon 
mass is the reduced mass, and <x is the Pauli matrix for 
the nucleon's spin. The function V ws is defined by 



V ws (r; V q ,k, R ,a,f3) 
N - Z 



-V 



1±K- 



A 



(G) 



1 + exp[dist s (r,/3)/a] 



where dists(r,/3) is the perpendicular distance between 
the point r and the surface S (taken with the minus sign 
inside the nucleus). 

For axially symmetric nuclear shapes, the surface £ is 
defined by 



R{9;R ,/3) = R c v {p) 



where (3 = {(3\} and c v ([3) is a normalization factor to 
conserve the volume. 

As the standard parameter set which corresponds to 
the actual nuclear potential, we employ the universal pa- 
rameter set given in Table 1 of Cwiok et al. [lj| for the 
main results. Precisely speaking, we do not use the depth 
of the central potential as it is, but modify it for each 
nucleus in a way explained in Sec. Ill CI and Ref. [TTI j . 
We also use five other parameter sets, which we call 
Wahlborn [2(1 Rost [H[, Chepurnovpl , Wyss-1 [Hj], 
and Wyss-2 [2J] (see Table I of Ref. [25| ) > to confirm the 
independence of the conclusions from the choice of the 
standard parameter set. 

Similarly to the Nilsson case, two constants f a and fi s 
are introduced by which to multiply the standard values 



of the surface diffuseness a and the spin-orbit potential 
strength A so , respectively. By changing the combina- 
tion of the values of these multipliers (f a , fi s ) from the 
standard values (1, 1), we examine how the prolate-shape 
dominance depends on the combined effects of the surface 
diffuseness and the spin-orbit potential. 

It should be noted that the centrifugal potential en- 
hances the effect of diffuseness on orbits having large an- 
gular momentum by pushing their wave functions out 
strongly on the surface. Thus, the effect of changing a 
by several tens of percents is quite large despite the fact 
that the standard value a = 0.7 fm is much smaller than 
the nuclear radius ~ 10 1 fm. 

The multiplier fi s of the WS potential is approximately 
equivalent to fi s of the Nilsson potential because both 
of them are the multipliers to the spin-orbit potential. 
Indeed, as we will show in Sec. IIIII the same value of fi s 
seems to correspond to the same situation between the 
two kind of potentials as far as the ratio of prolate-shape 
nuclei concerns. 

On the other hand, f a of the WS potential and /;/ of 
the Nilsson potential work in different ways to change 
(literally or effectively) the radial profile of the poten- 
tial. A finite- height cavity (i.e., square- well) potential is 
obtained in the limit f a — > of the WS potential and it 
can also be approximated by the Nilsson potential with a 
certain value of (> 1, dependent on the mass number 
A). A harmonic oscillator potential (oc r 2 ) is obtained 
by setting = of the Nilsson potential and it can 
also be approximately expressed by a certain value of f a 
(> 1, depending on A) for the WS potential. We will 
show in Sec. IIII Al that f a = 0, 1, and 2 correspond very 
roughly to = 1.5, 1, and 0, respectively, as far as the 
ratio of prolate-shape nuclei concerns, which is an aver- 
age over Oxygen to superheavy elements. One may also 
think roughly that WS potentials with f a > 2 cannot be 
expressed by Nilsson potentials while Nilsson potentials 
with /;/ > 1.5 or fa < cannot be expressed by WS 
potentials. 

In this paper, we calculate the ratio of prolate nuclei for 
31 x 19 = 589 kinds of artificially modified WS potentials 
defined by the combination of values of the multipliers 
(fis,. f a ) in ranges -1.5 < fu < 1.5 and 0.2 < f a < 2 
with Af ls = Af a = 0.1. 

Finally, we remark on the dependence of spin-orbit 



splittings, e - 



on the surface diffuseness. 



Concerning the Nilsson potential, the effects of two mul- 
tipliers fi s and are orthogonal and fu never affects 
spin-orbit splittings. As for the WS potential, however, 
spin-orbit splittings can be affected by f a through the 
radial form factor W W s in V S o of Eq. ©. (All the WS 
parameter sets postulate a so = a CE .) However, the net 
effect turns out to be small on the average after the in- 
tegral is done over the radius. 

In Fig.[IJ we show the average value of spin-orbit split- 
tings in 208 Pb as a function of f a . The potential is spher- 
ical and the parameter set is the universal set. In cal- 
culating the average, we consider all the pairs of levels 
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FIG. 1 : Average spin-orbit splitting in units of MeV of 208 Pb 
versus a multiplier f a to the surface diffuseness for six sets of 
parameters of the WS potential. 



The universal parameter set is exceptional in the sense 
that R so /R C e < 1 for neutron and > 1 for proton. We 
show the average over neutron levels and that over pro- 
ton levels separately for the splittings calculated with 
this parameter set in Fig. [5J The curve labeled as "neu- 
tron+proton" is the same as the curve for this parame- 
ter set in Fig. [T] One can see that the neutron splitting 
decreases by 25% while the proton splitting increases by 
100% in the plotted interval of f a - It is not clear whether 
the effects cancel between neutrons' decreasing and and 
protons' increasing splittings concerning the quantities 
like the ratio of prolate nuclei. From this point of view, 
the universal parameter set may not be the most suit- 
able set for the purpose of this paper. However, it is not 
a major demerit since the principal results of this paper 
like Fig. [3] and Fig. [T5] do not look very different between 
any parameter sets. Anyway, we regard the universal pa- 
rameter set as the standard and mainly show the results 
calculated with that parameter set. 



whose energies are negative throughout the plotted re- 
gion 0.1 < f a < 2. The number of states in each pair 
of levels (4/ + 2) is used as the weight in the averaging 
procedure. We do not distinguish proton and neutron 
levels for this figure. For the Wahlborn, Chepurnov, uni- 
versal, and Wyss-1 parameter sets, the variation in the 
average splitting is small (~ 10%). For the Rost and 
Wyss-2 parameter sets, the variation is larger (monoton- 
ically decreasing). We find that these different behaviors 
are originated in the difference between R so and R C e, i-e., 
the splitting does not change very much if R so = R CE , 
while it is a monotonically decreasing (increasing) func- 
tion of f a if i? so < R CE (R so > Rce)- The spin-orbit 
splitting of each pair of levels (not shown) behaves in the 
same way as the average (shown) except when the upper 
level is close to zero energy. 
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FIG. 2: Spin-orbit splittings in units of MeV of 208 Pb av- 
eraged over neutron levels, proton levels, and all the levels 
versus f a calculated with the universal parameter set for the 
WS potential. 



C. The Strutinsky method and its extensions 

The shape and energy of the ground state of each nu- 
cleus are determined by the shell correction method [|| . 
However, the conventional method cannot be applied to 
finite-depth potentials owing to the continuum part of 
the spectrum [2(J. In order to cope with this difficulty, 
we employ an improved shell correction method which we 
have recently proposed (the reference density method Q) 
based on the Kruppa prescription . 

In shell-correction approaches, the total energy of a nu- 
cleus is assumed to be expressed as a sum of the macro- 
scopic and microscopic parts, 

E = E mac + E mic = -Eldm + (Ebcs ~ -^bcs), (8) 

where .E m ac = -Eldm is the energy of the liquid- drop 
model, whose parameters are taken from Ref. [27|. The 
energies, Ebcs and -Ebcs, are calculated as 




where the occupation probability v 2 (e) and the pairing 
gap A (v 2 (e) and A) are obtained by solving the BCS 
equation for the seniority-type pairing interaction, whose 
matrix elements are a constant G, with discrete level den- 
sity g(e) (Strutinsky smoothed level density 5(e)). 

The interaction strength G is determined such that the 
smoothed pairing gap agrees with the empirical smooth 
trend, 

A = / A x 13 A- 1 ' 2 (MeV), (10) 

where the multiplier /a is one (which corresponds to 
actual nuclei) for most of the calculations except in 
Sec. IIIIBl where we try three other values ( /a = 0. 0.7, 
1.2). 
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The depth of the central potential V GE is adjusted for 
the consistency between the microscopic (in the Thomas- 
Fermi approximation) and macroscopic Fermi energies 
for spherical shapes neglecting the spin-orbit potential. 
This adjustment is necessary in order to treat nuclei far 
from the /3 stability line because the available parameter 
sets of the WS potential do not reproduce the location 
of the drip lines predicted by the macroscopic part of 
energy. This adjustment is also indispensable when f a ^ 
1, because the Fermi level is strongly affected by f a . 

More details on the methods of calculations are de- 
scribed in Ref. @j. 



D. Numerical precision 

We have confirmed that numerically calculated single- 
particle spectrum of the WS potentials is sufficiently pre- 
cise to calculate the ratio of prolate nuclei. It should be 
noted that usually recommended size and frequency of 
the oscillator basis are not guaranteed to be sufficient to 
diagonalize Hamiltonians when potentials arc artificially 
modified, i.e., when fi s or f a is not one. 

We prepare the basis of diagonalization as cigenfunc- 
tions of a Hamiltonian of a nuclcon of mass m in a har- 
monic oscillator potential, whose frequency is determined 
in such a way that the major shell spacing is equal to 
Ufa = fu X 41 A- 1 / 3 MeV. Here, f u is an adjustment 
factor, usually takes on 1.2. When f a > 1, the spatial 
extension becomes larger so that the oscillator length of 
the basis oc f u should be increased, i.e., / w should be 
decreased. We determine empirically f u as a function of 

fa j 

f 1.2 (0 < /„ < 0.5), 

f u = { 1.4-0.4/ o (0.5 < f a < 1.5), (11) 

{ 0.8 (/„ > 1.5), 

so as to obtain more precise results with smaller basis. 

The single-particle basis is truncated according to the 
maximum number of the oscillator quanta N osc . In this 
paper, we diagonalize mostly in a subspace N osc < 20. 
Concerning the calculations for Figs. [31 HI QTJ and [131 
we consider only iV osc < 16 because of the required large 
computations. As we show in Sec. IIII A[ the results are 
quite close to each other between the calculations with 
N osc < 16 and N osc < 24 as far as the ratio of prolate 
nuclei concerns. 



E. Set up of the numerical calculations 

For each potential specified by a set of multipliers 
(fis,fu) or (fi s , f a ), we calculate the ground-state shape 
for all the even-even nuclei with 8 < Z < 126 and 
8 < N < 184 between the proton and neutron drip 
lines predicted by the Bethe-Wcizackcr (macroscopic) 
mass formula (2148 nuclei). To find the ground-state 



shape for each of these 2148 nuclei with each of the 
589 kinds of the WS potentials, we calculate the nuclear 
total energy surface versus the quadrupole and hexade- 
capolc deformation parameters ($2, Pa) at 51x21 = 1071 
points in ranges —0.5 < /?2 _ 0.5 with A/?2 = 0.02 and 
-0.3 < /3 4 < 0.3 with A/3 4 = 0.03. For the Nilsson- 
Strutinsky calculations of the 496 kinds of the Nilsson 
potentials, (3\ should be read as e\. 

The scale of required numerical computation is quite 
large. In order to obtain each of the six panels Fig. [3] (a) 
and Fig. [T5] (a-e) corresponding to the different param- 
eter sets of the WS potentials, we have to calculate the 
total energy 589 x 2148 x 1071 = 1.4 x 10 9 times, each 
time of which is for a nucleus using a given-shape po- 
tential. We have performed such computations in 64 
threads in a PC cluster having eight Intel Core i7-920 
CPU (2.66GHz, four physical = eight logical cores in a 
CPU). It has typically taken 28 hours for a physical core 
to complete a nuclear chart for a given combination of 
(fis,fa), and 20 days for the cluster to complete all the 
combinations, when we treat the WS potential in an os- 
cillator space truncated by N osc < 16. Compared with 
the WS potentials, it takes negligibly small time to per- 
form the corresponding Nilsson-Strutinsky calculations. 



F. Definition of the ratio of prolate nuclei 

After completing all the Strutinsky-method calcula- 
tions, we count the number of prolate, oblate, and spher- 
ical nuclei in each nuclear chart to calculate the ratio of 
prolate nuclei R p as a measure of the prolate-shape dom- 
inance. We use two different definitions of the ratio in 
order to demonstrate that the main conclusions of this 
paper are not affected by the choice of the definition. 

By N , N s , and N p , let us denote the number of nuclei 
whose ground states are oblate (fa < —0.05), spheri- 
cal (Ifol < 0.05), and prolate (fo > 0.05), respectively, 
among the 2148 even-even nuclei defined in Sec. Ill El with 
those having positive Fermi levels excluded. Using these 
numbers, the ratio of prolate nuclei may be defined as 



In Ref. |5|, however, we have proposed a different defi- 
nition in order to exclude shape transitional nuclei which 
often have an almost flat bottom in the nuclear energy 
curve ranging typically between — 0.2 < fa < 0.2. It is 
better to exclude them because precise ground states of 
such nuclei are superpositions of prolate and oblate shape 
states owing to large zero-point quantum fluctuations. 
There is little point whether the minimum happens to be 
located in the prolate or oblate side. After an extensive 
examination of the landscapes of total energy curves of 
all the even-even nuclei, we found that excluding spheri- 
cal and transitional nuclei is nearly equivalent to includ- 
ing only those nuclei which have both prolate and oblate 
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minima (discarding very shallow ones [5[)- Hence we have 
defined alternatively the ratio of prolate nuclei as 



(a)WS 



Rn 



p 



P 



(13) 



where N™ ( N° 



) denotes the number of nuclei which 
has both oblate and prolate minima and the prolate 
(oblate) minima has a lower energy. (Superscript "DM" 
stands for "double- minimum" ) . 

In Sec. IIII A[ we will demonstrate that the main con- 
clusions of this paper is not altered essentially whether 
we use either R p or R' p . The advantage of R p is the ex- 
clusion of transitional nuclei while the merit of R' p is only 
simplicity. Hence we mainly use R p in this paper. 



III. RESULTS 

Ratio of prolate nuclei versus potential 
parameters 
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FIG. 3: The ratio of prolate nuclei R p for the WS (a) and the 
Nilsson (b) potentials. The abscissa is the multiplier /; s to 
the spin-orbit potential for both panels. The ordinate is the 
multiplier f a to the surface diffuseness a of the WS potential 

(a) and the multiplier fu to the I 2 potential of the Nilsson 
potential (b). The universal parameter set is used as the 
standard one for the WS potential. 

In Fig. 02 the ratio of prolate nuclei R p defined by 
Eq. (jTB")) is shown as a function of fi s and f a for the WS 
potential (a) and /;/ and fi s for the Nilsson potential 

(b) . We use the universal parameter set [lj| for the WS 
potential. Fourth order Legendre polynomials in fi s and 



FIG. 4: Same as Fig. [3] but for a different definition of the 
ratio of prolate nuclei R' p . 



those in f a (or fu) are used for the interpolations to draw 
the contours. 

In Fig. 21 the ratio of prolate nuclei R' p of a simpler 
definition of Eq. (|12[) is shown in the same manner as 
R p in Fig. [3l One can see the behavior of R' p is quite 
similar to that of R p . It can also be confirmed in detail 
in Figs. [SHU which show curves for R p and R' p in various 
cross sections of the panels in Figs. [3] and |U The distance 
between the two curves is mostly ~ 0.05. The largest 
discrepancy is 0.11 at (fu, fu) = (0,1) in Fig. Mjo)- We 
have also checked that a third definition 



R „ = ^p 
p N p + N s + N 



(14) 



does not change the behavior very much, either (no fig- 
ures). Therefore, the following discussions are indepen- 
dent of the definition of the ratio of prolate nuclei. 

Result for the Nilsson potential shown in Fig. [3Jb) is 
essentially the same as Fig. 1 of Ref. Q , except that the 
area of fu < is not shown in this paper. 

In a line = (Fig. EJa)), i.e., when there is no 
spin-orbit potential, R p (R' p ) increases monotonically as 
a function of /;/, from 60% (56%) at fu = (i.e., the 
harmonic oscillator potential) to 80% (71%) at /;/ = 1 
(corresponding to the actual surface thickness). This be- 
havior agrees with Frisk's argument [|[ that a cavity(- 
like) potential prefers prolate shapes more than oblate 
ones. 

In a line fu = 1 (Fig. EJb)), i.e., with the surface 
diffuseness of actual nuclei, R p as well as R p oscillates 
strongly as a function of fi s . It takes local maximum val- 
ues at fi s ~ ±1 and local minimum values at fi s ~ ±0.5. 
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FIG. 5: The ratio of prolate nuclei R p and R' p for the WS 
potential for (a) f a = 0.4, (b) f a = 1, and (c) f a — 1.6 versus 
the multiplier to the spin-orbit strength fi s . 
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FIG. 6: The ratio of prolate nuclei R p and R' p for the Nilsson 
potential for (a) fu = 1.5, (b) fu = 1, and (c) fu = versus 
the multiplier to the spin-orbit strength / is . 



It means a strong interference between the effects of the 
Is and l 2 potentials. Minima of R p are located close to 
this line fu = 1 where R p ~ 40%, i.e., oblate shapes 
dominate over prolate shapes if the spin-orbit potential 
is weakened by 50% while keeping the surface diffuseness 
of the potentials at the standard value. 

The result for the WS potential is shown in Fig. EHa). 
[1H It looks quite similar to Fig. EJb) for the Nilsson 
potential. 

In a line fi s = (Fig. Ufa)), R p (R' p ) increases almost 
monotonically as a function of —f a> from 68% (64%) at 



f a = 2 to 91% (86%) at f a = 0.2. This behavior during 
the transition from the harmonic oscillator to a cavity- 
like potentials agrees qualitatively with that of the Nils- 
son potential and with Frisk's argument Q. Quantita- 
tively, however, R p = R' = 73% at f a = 1 (actual nu- 
clear central potential) is not sufficiently dominant. The 
situation is worse in a line fi s = 0.5 (Fig. E£b)). Only by 
the assist of the actual-strength spin-orbit potential, i.e., 
fis = 1 (Fig. [7£c)), R p can exceed 80% when f a = l. 

As for the Nilsson potential, this assist of the spin-orbit 
potential works stronger to change the landscape of the 
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FIG. 7: The ratio of prolate nuclei R p and Tip for the WS 
potential for (a) fu = 0, (b) fu = 0.5, and (c) fu — 1 versus 
the multiplier to the surface diffuseness f a . 
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FIG. 8: The ratio of prolate nuclei R p and R' p for the Nilsson 
potential for (a) fu = 0, (b) fu = 0.5, and (c) fu = 1 versus 
the multiplier to the I 2 term fu. 



curve. When fi s = 0.5, it generates a large- amplitude 
oscillation versus fu (Fig.[8]Jb)), whose minimum at fu = 
0.95 corresponds to the minimum at (fi s ,fu) = (0.45, 1) 
in Fig.[6jb). It is this oscillation that has moved two local 
minima upward in two-parameter plots of R p and of R' 
compared with those for the WS potential: For the WS 
potential (Fig. [3] (a)), the minima are located around a 
line f a ~ 1.6, not f a = l, while for the Nilsson potential 
(Fig.^b)), they are just in a line fu = 1. The amplitude 
of the oscillation in a line fu = 1 (Fig. [H^b)) is larger 
than that in a line f a = 1.6 (Fig. Efc)) or that in a line 



f a = 1 (Fig. E^b)) by a factor ~ 3. The reason for the 
exaggeration of the amplitude of the oscillation by the 
Nilsson potential may be the approximate treatment of 
the radial profile of the potential in terms of the I 2 term 
and the affinity between this I 2 term and the Is term in a 
sense that both of them contain the same orbital angular 
momentum operator I. Anyway, we have confirmed a fact 
that the interference also exists for the WS potential. 

In Fig. [SI we show that our single-particle basis is suffi- 
ciently large. Panel (a) is the same as Fig.[5][b) and panel 
(b) is the same as Fig. [7f a) except that results with three 
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FIG. 9: Dependence of the ratio of prolate nuclei R p on the 
size of the oscillator basis specified by N osc , which is the max- 
imum oscillator quantum number of the single-particle basis. 
In panel (a), R p is plotted versus the multiplier to the spin- 
orbit potential fi„ while the surface diffuseness is frozen at the 
standard value (f a = 1). In panel (b), R p is plotted versus 
the multiplier to the surface diffuseness f a while the spin-orbit 
potential is kept turned off (/;„ = 0). 
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FIG. 10: The dependence of the ratio of prolate nuclei R p 
on the pairing strength for the WS potential. Four kinds of 
curves correspond to different values of a modification factor 
/a which determines the strength of the paring force. For 
/a = 1, the empirical trend is used. For /a = 0, there is no 
pairing correlation. Panel (a): as a function of the spin-orbit 
strength /; s with a normal diffuseness f a = 1. Panel (b): as a 
function of the diffuseness f a without a spin-orbit potential. 



different sizes of basis are compared. One can see in both 
panels that the curve for N osc < 16 is quite close to the 
curve for iV osc < 24. 



B. The effect of the pairing strength 

In Fig. [TUl the ratio of prolate nuclei R p is shown for 
four values of the multiplier /a appearing in Eq. (|10p for 
the WS potential. The pairing strength can be controlled 
with this multiplier. 

In panel (a), R p is shown as a function of /; s while 
the surface diffuseness is the standard one (f a = 1). The 
short dash curve labeled as /a = 1.0 is the same as the 
solid curve in Fig.[5jb). One can see that a stronger pair- 
ing increases the amplitude of oscillation. The same ten- 
dency has also been found for the Nilsson potential [111 ]. 
The change is largest at prolate-shape dominant points 
°f fls ~ ±1, where R p is increased further from a high 
value to a higher value. This can be understood that gen- 
erally shallower minima in oblate nuclei are more easily 



dissolved and merged with the other minima to form a 
single spherical minimum than generally deeper minima 
in prolate nuclei. 

In panel (b) , R p is drawn versus f a while the spin-orbit 
potential is kept turned off (/; s = 0). The short dash 
curve labeled as /a = 1.0 is the same as the solid curve 
in Fig. Eta). One can see an overall trend that stronger 
pairing correlation increases the magnitude of the slope. 
Pairing correlation increases R p at / a < 1, which may be 
understood by the same interpretation given in the last 
paragraph. Pairing correlation decreases R p at 1 < f a < 
1.5. At f a > 1.5, curves are reordered and the trend of 
monotonic increase of R p versus — f a is reversed for /a > 
1. These behaviors require a different interpretation yet 
to be explored. 

C. Mass dependence of the ratio of prolate nuclei 

We decompose Fig. G^a) into three mass regions, A < 
100, 100 < A < 200, and A > 200 to obtain the three 
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FIG. 11: Same as in Fig. [3ja) but decomposed into a light 
A < 100 (a), a medium 100 < A < 200 (b), and a heavy 
200 < A (c) mass regions. 



panels of Fig. [TT] One can see that the patterns in these 
three mass regions are rather different; the numbers of 
oscillations in both directions (/ Q and /; s ) are larger in 
the heavier mass region. Namely, the ratio of prolate 
nuclei changes more rapidly when the potential is modi- 
fied. This pattern change can be understood by the fact 
that the level density is larger in heavier system, where 
the prolate- and oblate-driving orbits coexist in a narrow 
energy range at the Fermi surface and a modification of 
the potential more easily changes the equilibrium defor- 
mation. Alternatively, one can think that the pattern 
is "shrunken" in two axes from the light to the heavy 
regions, although new pronounced low R p areas appear 
in the range of small f a < 0.5 in the heavy mass region 
(Fig.QjJc)). The shrinkage in the f a axis may be ascribed 
to the ratio of a constant thickness a to the nuclear ra- 
dius oc A 1 / 3 , while that in the /; s axis may be related to 
the ratio of the spin-orbit splitting oc A osc oc A 1 / 3 to the 
major shell spacing oc A^ 1 ^ . This shrinkage may cause 
a blurring of the pattern averaged over all the mass re- 



gions. It may be possible that, by employing a properly 
A-scaled parametrizations for the surface diffuseness a 
and the spin-orbit potential strength A so , instead of sim- 
ple constant multipliers f a and fi s , one obtains a little 
sharper pattern than Fig. [3J 



D. Absolute number of nuclei of prolate, oblate, 
and spherical shapes 

In Fig. [T^J we show not ratios but absolute numbers 
of nuclei of prolate, oblate, and spherical shapes at their 
ground states, i.e., N p , N Q , and N s defined in Sec. Ill Fl 
The numbers of nuclei are plotted versus fi s while f a is 
frozen at 1. Because these numbers depend on the mass 
number more strongly than ratios, we discuss separately 
the light, medium, and heavy mass regions introduced 
in Sec. IIII CI As it is discussed in the previous subsec- 
tion, one observes considerable differences in three mass 
regions; the three numbers, 2V P , N Q , and N s , oscillate 
more frequently in the heavy mass region. In Fig. [T2l a') , 
which shows the numbers of nuclei in the light mass re- 
gion, the oscillation in N p for |/; s | < 0.8 can be ascribed 
to shape transitions of the ground state of each nuclei be- 
tween prolate and spherical shapes because N a is almost 
constant in that interval of fi s . In Fig. Il"27 b) , which is 
for the medium mass region, the amplitude of oscillation 
in N p is most pronounced among the three mass regions. 
The oscillation is attributed evenly to shape transitions 
between prolate and spherical shapes and those between 
prolate and oblate shapes. In Fig. H2Ic), which corre- 
sponds to the heavy mass region, the amplitude of the 
oscillation in N p is reduced. This may be a consequence 
of the fact that, from the light to the heavy regions, the 
prolate nuclei continues to become more and more dom- 
inant over spherical and oblate nuclei for the normal dif- 
fuseness f a = 1. 



E. Map of deformation on the nuclear chart 

In Fig. 1131 we plot the quadrupolc deformation pa- 
rameter /?2 of the ground states of even-even nuclei in 
the (N, Z) plane. In panels (a-c), the WS potential is 
used. The parameter sets are the universal parameter 
set which is the choice of this paper (a), the Chepurnov 
parameter s et |22l| which is one of the classical (b), and 
the Wyss-2 [24[ parameter set which is one of the latest 
(c). In panel (d), the Nilsson potential is used and e 2 is 
plotted instead of /?2- Comparing the predictions of these 
potential parameter sets, one sees that the locations of 
prolate and oblate nuclei are roughly the same among 
the four panels. This means that the conclusions of this 
paper is not restricted to the mainly employed universal 
parameter set. 

With a more careful look at these maps, one finds that 
a rectangular area 28 < Z < 40, 50 < N < 82 is the 
area where the largest number of nuclei have differently 
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FIG. 12: The number of prolate (JV P ), oblate (JV ), and 
spherical nuclei (N a ) in a light (a), a medium (b), and a heavy 
(c) mass regions as a function of /;„ (with /„ = 1). 



in the A ~ 80 region, see e.g. Ref. [28|]. Compared to 
the A ~ 80 region, this rectangular area is much wider 
so that it may be a more suitable region to study the 
issue of the shape coexistence and/or the shape phase 
transition in the light of recent developments of unsta- 
ble beam facilities. It should be noticed that a smaller 
area 28 < Z, N < 40 has a similar characteristic, which 
provides t ypic al examples of the shape coexistence phe- 
nomenon [291 ] . where the spherical shape also play an 
important role. 

The deformation map changes drastically when the 
multipliers are changed. In Fig. [H] we plot the 
quadrupole deformation predicted with the universal 
parameter set for the WS potential modified with a mul- 
tiplier f u = 0.5 (a), (b) ,-0.5 (c), and -1 (d) while 
f a is frozen at 1. Comparing these four panels together 
with Fig. HUT a) for fi s = 1 , one sees a substantial change 
of the distribution of the shapes caused by varying the 
value of fi s . 

One also notices that positive and negative values of 
fls of the same magnitude (i.e., Fig.[T3la) and Fig. [14T d). 
Fig. IT47 a) and Fig. H4T c)) lead to very similar distribu- 
tions. 

For fi s — ±1, oblate-shape nuclei are located in a nar- 
row region just below major-shell closures. On the other 
hand, for /; s = ±0.5, they occupy a large area in the lat- 
ter half of a major-shell filling. The largest of such areas 
are in 50 < Z < 80, 90 < N < 120 in Fig. [Tg^a.c). 

It is interesting that parameters fi s = and /; s = ±1, 
which lead to similar levels of prolate-shape dominance, 
have completely different shape distribution maps. This 
seems mainly because, when the spin-orbit potential is 
missing ( fi s — 0), subshell closures are strong enough to 
subdivide conventional major shells. 



F. Dependence on WS parameters 

In Fig. [151 we show the ratio of prolate nuclei R p for 
various WS parameter sets. Combined with Fig. [3] (a), 
all the six sets of the WS potentials are investigated. In 
all the figures the basic pattern of the contour plot is 
the same, although the precise positions of the low R p 
regions are slightly shifted in each case. Especially, the 
physical points, fi s = f a = 1, is near one of the highest 
points, indicating the strong prolate-shape dominance for 
all parameter sets. By comparing them, one can say that 
the conclusions of this paper are not essentially altered 
by changing the standard potential. 



predicted shapes by different parameter sets; the predic- 
tion of the Nilsson potential is also very different from 
those of the WS potentials. In this nuclei of large 

prolate deformations and those of large oblate deforma- 
tions are located side-by-side directly, meaning the first 
order phase transition in the shape. A similar situation 
that the distinct prolate and oblate minima coexist in 
the same energy region has been extensively discussed 



G. Single-particle level densities 

Wc give an example to illustrate the way how the po- 
tential parameter changes the shape of a nucleus through 
the change of the level density. We choose a nucleus 
1 ggErio4, which is located near the center of a region 
52 < Z < 80, 90 < N < 110, where the ground state 
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FIG. 13: Quadrupole deformation parameter of the ground 
states of even-even nuclei calculated with standard parame- 
ters for the WS potential (a-c) and the Nilsson potential (d). 
The used parameter set for the WS potential are the universal 
(a), the Chepurnov (b), and the Wyss-2 (c). 
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FIG. 14: Quadrupole deformation parameter of the ground 
states of even-even nuclei using the universal parameter set 
for the WS potential with multipliers f a = 1 and fi s — 0.5 
(a), (b), —0.5 (c), and —1 (d). Horizontal and vertical dot 
lines indicate the magic numbers of actual nuclei, not of the 
modified potentials. 



is prolate for the standard potential and oblate for a 
modified potential with (fi s ,f a ) = (0.5,1) as shown in 
Fig. [Ha). 

In Fig. [T|)l the oscillating part of the neutron's level 
density defined by, 

dg(e) = <? 7 =o.5(e) - <? 7 =i.2(e), (15) 
of the ground state of this nucleus is plotted as a func- 



tion of the multiplier to the spin-orbit potential /; s (ab- 
scissa) and the single particle energy in unit of fko (or- 
dinate). The Fermi level for neutrons is designated with 
a black curve. The quantity 7 in Eq. (|15[) is the Struti- 
nsky smoothing parameter in unit of fwj, and the value 
7 = 0.5 is chosen to clearly show the oscillation due to 
the shell effect relative to the average level density with 
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FIG. 15: Ratio of prolate nuclei R p calculated by using dif- 
ferent WS parameter sets as the standard, the Wahlborn (a), 
the Rost (b), the Chepurnov (c), Wyss-1 (d), and Wyss-2 (e). 



7 = 1.2, with which all Strutinsky smoothed quantities 
are calculated. 



Fig. TOa) shows the level density at a prolate shape 
{Pi = 0.282), in which the level density at the Fermi level 
is low (high) at fi s = 1 (0.5). Fig. [TfJlb') shows the level 
density at an oblate shape (/?2 = —0.220), in which the 
level density at the Fermi level is high (low) at fi s = 1 
(0.5). The change from a prolate to an oblate shape of 
this nucleus due to the change of fi s from 1 to 0.5 can 
be explained in this way in the macroscopic-microscopic 
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FIG. 16: Oscillating part of the neutron's single-particle level 
density 8g{e) of 172 Er as a function of the multiplier to the 
spin-orbit potential fi s and the single-particle energy divided 
by hw. The shape of the potential is prolate in panel (a) and 
oblate in panel (b). The Fermi level is designated with a black 
curve. 



theory. 



H. 



Prolate, oblate, and spherical ratios as 
functions of Fermi energy 



As it is explained in detail in Sec. Ill El we have calcu- 
lated the ground state masses of 2148 even-even nuclei in 
the nuclear chart for each modified set of parameters of 
the WS potential. In view of the results, there are many 
interesting findings apart from the prolate-shape domi- 
nance. Among them, we present the result on how the 
deformation of unstable nuclei changes compared with 
the stable nuclei. In Fig. [T71 three ratios of the prolate, 
oblate, and spherical nuclei, i.e., Rp defined by Eq. (p~4|) 
and 



R'' = 



N 



N v + N s + N 



R'J 



iVs 



N D + N S + N 



(16) 



are shown as functions of the neutron (panel (a)) or pro- 
ton (panel (b)) Fermi energy. Note that these ratios are 
now defined with including N s in contrast to the defini- 
tions of R p and R p used in previous subsections which 
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FIG. 17: Ratios of prolate, oblate, and spherical nuclei as 
functions of the Strutinsky smoothed Fermi energy A for the 
neutrons (a) or for the protons (b). The parameter set is the 
standard (/; s = f a = 1) universal set. 



exclude it. For every one MeV interval of the calculated 
Fermi energy, we count the numbers of nuclei N p , N a , 
and TVs, calculate their ratios R p \ i?", and i?", and plot 
the ratios with the abscissa at the center of the inter- 
val in Fig. 1171 The Strutinsky smoothed Fermi energy 
is used as the calculated Fermi energy. There are a few 
other choices for the Fermi energy; for examples, half of 
the two nucleon separation energy, £211/2 (5^/2), or the 
chemical potential obtained by the BCS calculation. The 
latter is not well-defined in the case of vanishing pairing 
gap. The former choice has been checked that the resul- 
tant plot is very similar to Fig. [T7] 

It has been sometimes discussed that in the weakly- 
bound system the wave functions of the orbits near the 
Fermi surface spread out and then the shell effects play a 
minor role. Hence, the ratio of spherical nuclei increases 
in weakly-bound unstable nuclei. Such a trend is clearly 
seen in Fig. ITTf a); the ratio of prolate nuclei decreases 
from 0.6 at A n w —8.5 to 0.4 at A n « —0.5 in balance 
with the increase of the spherical ratio. Even the prolate- 



shape dominance is inverted in the limit of weak-binding 
for neutrons, — A n — > 0. As for protons, there is no such 
trend; the ratios are almost constants as the Fermi energy 
is increased and the spherical ratio of nuclei even slightly 
decreases. This is because the existence of the Coulomb 
barrier prevents the weak-binding situation even in the 
limit of vanishing Fermi energy. 



IV. SUMMARY 

We search for the origin of the prolate-shape domi- 
nance of nuclear ground state deformation in (the combi- 
nations of) the properties of the single-particle potential 
for nucleons. We employ the Woods-Saxon (WS) poten- 
tial and apply to it the macroscopic-microscopic theory 
modified for the treatment of the continuum part of the 
spectrum. We change the surface thickness, the strength 
of the spin-orbit potential, and the strength of pairing 
correlations to examine their influences on the ratio of 
prolate nuclei over more than two thousand even-even 
nuclei. We observe strong interference between the ef- 
fects of the surface thickness and the spin-orbit poten- 
tial, an especially interesting consequence of which is an 
oscillation versus the spin-orbit potential strength. We 
also find that pairing correlations enhance prolate-shape 
dominance. 

These results are compared with the results of the Nils- 
son potential to elucidate some of the special features of 
the Nilsson potential, especially an exaggeration of the 
amplitude of this oscillation for the ratio of prolate nu- 
clei. 

We repeat calculations for six different parameter sets 
of the WS potential. We find these parameter sets can 
be classified according to the relation between the radii 
of the central and the spin-orbit potentials, which deter- 
mines the effect of the surface diffuseness on spin-orbit 
splittings. However, as far as the behavior of the ratio of 
prolate nuclei concerns, different parameter sets do not 
change the conclusions of this paper substantially. 

The difference of the definition of the ratio of prolate 
nuclei is also shown to be unimportant. The error in 
the ratio of prolate nuclei due to the truncation of the 
oscillator basis to express single-particle wavefunctions 
is shown to be sufficiently small. 

Maps of ground state quadrupolc deformation arc 
shown and discussed in detail for the physics behind the 
ratio of prolate nuclei, which is only a single number. 
An example is given to illustrate the relation between 
the level density and the shape of the ground state in 
the macroscopic-microscopic theory. The ratios of pro- 
late, oblate, and spherical nuclei arc investigated also as 
functions of the Fermi energy. It is found that the spher- 
ical ratio increases and even exceeds the prolate ratio in 
the limit of weak-binding for the neutron rich unstable 
nuclei. 
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